The Fractional Quantum Hall Effect of Tachyons in a Topological Insulator Junction 
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We have studied the tachyonic excitations in the junction of two topological insulators in the pres- 
ence of an external magnetic field. The Landau levels, evaluated from an effective two-dimensional 
model for tachyons, and from the junction states of two topological insulators, show some unique 
properties not seen in conventional electrons systems or in graphene. The v = | fractional quantum 
Hall effect has also a strong presence in the tachyon system. 
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The surface state of recently discovered three- 
dimensional topological insulators Q contains a single 
Dirac cone and as a result, the charge carriers on the sur- 
face are characterized as massless Dirac fermions. Some 
of the properties of these particles are well known from 
another topological system, the graphene 3]. We have 
shown earlier [J] that the dispersion relation of the sur- 
face excitations in a junction between two such topo- 
logical insulators (TIs) show several very unique prop- 
erties. Most importantly, under certain conditions these 
excitations exhibit tachyon-like dispersion relation 0, Q 
corresponding to superluminal propagation of the Dirac 
fermions along the interface of the two TIs. In addition 
to the tachyonic dispersion, the junction states can also 
support the usual massless relativistic dispersion relation 
of the Dirac fermions Q • Here we report on the proper- 
ties of these tachyonic excitations in an external magnetic 
field. We discuss the unique nature of the Landau levels 
(LLs) of these tachyons and the interaction properties of 
tachyons in a strong magnetic field, which leads to the 
fractional quantum Hall effect (FQHE) of tachyons. 

Effective two-dimensional (2D) model of tachyonic ex- 
citations: The tachyons in the junction of two TIs can be 
described by the effective 2D matrix Hamiltonian of the 
Dirac fermions but with imaginary Fermi velocity [H (in- 
stead of the imaginary proper mass commonly attributed 
to the tachyons @|, e.g., to the neutrinos Q) 
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where p± = p x ± p y is the 2D momentum, A is 
the effective mass of the tachyons, i«j is the imaginary 
Fermi velocity, and a are the Pauli matrices, corre- 
sponding to the spin degree of freedom of an electron 
(tachyon). The effective Hamiltonian thus described has 

the tachyon-like dispersion, E Tach (p) — ±-^/Aq — v Y p 2 , 

i 

where p = (p 2 . + p 2 ^ 2 . Typical values of the parameters 
are A ~ 0.3 eV and v l ~ 10 6 m/s. 

We subject the tachyonic system to an external mag- 
netic field, B, pointing along the z-direction, i.e., perpen- 
dicular to the junction between the TIs. The Hamilto- 
nian describing such a system in a magnetic field can 
be found from the Hamiltonian ((T|) by replacing the 



tachyonic momentum, (p xl p y ) by the generalized mo- 
mentum, (TT x ,TT y ) and introducing the Zeeman energy, 
A Z (B) = ^g s fi B B. Here fi B is the Bohr magneton and 
g s ~ 8 is the effective (/-factor of the tachyonic excita- 
tions, which we assume to be equal to the g-factor of the 
TI surface states 0, |T(| ■ The effective tachyonic Hamil- 
tonian in a magnetic field then takes the form 



H Tach (B) = [A + A Z (B)} a z + iv^an). 
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The LL energy spectrum corresponding to the Hamilto- 
nian [Eq. is characterized by an integer number n > 
(the LL index), and is defined as 



E n=Q = A Q + A Z (B) 
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Here s = ±1, £ = y 'eh/cB is the magnetic length, 
and we introduced the effective magnetic field B*(B) = 

^-{(A + A z (B))/v 1 ] 2 . 

The wave functions corresponding to the LLs [Eq. 
can be expressed in terms of the functions <p n m , which 
are wave functions of the conventional ('non-relativistic') 
LLs with index n and the intra-Landau level index, to, 
for example, the z component of the angular momentum. 
The LL wave functions of the tachyons are 



T (Tach) _ 



y 0,m 





(Tach) 
ro> l,s 



cos a 



n.s rn,m 



n,s Tn—l,m 



sy/1- nB/B' 



1/2 



(4) 



where a ns = sin 

The LL energy spectrum obtained from Eq. (j3j) is 
shown in Fig. 1. The energy spectrum for tachyons ex- 
hibits a few distinct features: (i) The LL energy spectrum 
is mainly restricted within the energy interval — A < 
E n s < A . (ii) At a given magnetic field B, only the LL 
with index n < B*/B can be observed. Therefore, for 
a given magnetic field, only a few LLs exist in the sys- 
tem and, for a large enough magnetic field (e.g., B > 32 
Tesla in Fig. 1(a)), only the n = LL survives. This 
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FIG. 1: (a) The LLs of the effective 2D tachyonic Hamilto- 
nian {TJ are shown as a function of the magnetic field. The 
numbers next to the lines are the LL indices. The n — 1 
LL wave functions at points A and B are identical to conven- 
tional 'nonrelativistic' LLs with indices and I, respectively. 
At point G, the LL wave function of the tachyon system is 
identical to the n = 1 LL wave function of graphene. (b) 
v — |-FQHE gap is shown for the n — and n = I LLs of 
tachyons. The FQHE gap is calculated for a finite-size system 
with N = 9 particles in a spherical geometry with parameter 
S = 12. The red and blue lines correspond to the n — 1 LL 
branches shown in panel (a). The FQHE gap is measured in 
units of the Coulomb energy, e c = e 2 /k£ q . The dashed curve 
is explained in the text. 



behavior of the LL of tachyonic excitations is totally dif- 
ferent from that in the 'non-relativistic' semiconductor 
systems with the LL energy spectrum E n oc nB, and 
"relativistic" graphene system with the LL energy spec- 
trum E n oc V nB Q. However, there is one similarity be- 
tween the tachyonic LL dispersion relations and that in 
graphene. In both cases there is one n = LL, the energy 
of which is independent of the strength of the magnetic 
field (without the Zeeman splitting). In graphene, the 
energy of this LL is E n=0 = 0, while for the tachyonic 
excitations, E n=0 = A . In both cases the corresponding 
wave functions are 4> n=0 m ■ 

To address the similarities and differences between the 
LL wave functions of the tachyonic system and those of 
graphene (or even conventional semiconductor systems), 
we consider the interaction properties of tachyonic ex- 
citations in a given LL. To characterize the strength of 
the inter-tachyonic interactions we study the strength of 
the FQHE, i.e., the magnitude of the FQHE gap. In the 
FQHE regime the electrons partially occupy a single LL 
and the properties of such a system are characterized by 
the inter-particle interactions within the corresponding 



LL [ll|] . The interaction strength within a given LL is 
determined from the Haldane pseudopotentials, V m [l2j], 
which are the interaction energies of two particles with 
relative angular momentum m. The pseudopotentials are 
determined from 12] 



p- q V(q)[F n (q)f L m {q*)e-<? , (5) 



where L m (x) are the Laguerre polinomials, V(q) — 
2ire 2 /(K,£ q) is the Coulomb interaction in the momen- 
tum space, k is the dielectric constant, and F n (q) is the 
form factor for the n-th Landau level. The form factor 
is determined by the structure of the LL wave functions. 
For the tachyonic system the form factor is given by 

F n=0 (q) = L {\q 2 ) (6) 
F n >i(q) = cos 2 a n s L n (\q 2 ) + sin 2 a ra>s L„_ 1 (\q 2 \l) 

The form factor of the n = LL [Eq. ©] is identical to 
that of graphene and also to that of the non-relativistic 
systems. However, for n > 1 the form factor of the tachy- 
onic system becomes unique. For graphene and for the 
non-relativistic systems the corresponding form factors 



are Fn 



(NR.) 



and Fi Gr) 



h(L n + respec- 



tively. In both cases the form factors are independent 
of the magnetic field. For the tachyonic system, on the 
other hand, the form factor (J7J depends on the mag- 
netic field, through the effective angle a n S (B). With in- 
creasing magnetic field the tachyonic form factor, F n>1 , 

changes from the non-relativistic value, Fn B:> [point B 
in Fig. 1(a)] or [point A in Fig. 1(a)], in a small 



magnetic field, B 

(Gr) 



0, to the form factor of graphene, 
F^> = \{L n + L n _ x ) [l]|, for B = B*. 

The FQHE with an incompressible ground state can 
be observed only in the LL with strong short-range re- 
pulsion, i.e., a fast decay of the corresponding pseudopo- 
tentials, V m , with m. Such a strong repulsion is realized 



only in the LL with a strong admixture of i 



'0,TO' 



There- 



fore, in a tachyonic system the FQHE is expected only in 
the n = and n = 1 LLs. To study the strength of the 
corresponding FQHE we numerically evaluate the energy 
spectrum of a finite A^-electron system in a spherical ge- 
ometry [12[ with the radius of the sphere VS£ . Here 
2S is the integer number of magnetic fluxes through the 
sphere in units of the flux quantum. For a given num- 
ber of electrons, the radius of the sphere determines the 
filling factor of the system. For example, the v = 1/q 
FQHE (q is an odd integer) corresponds to the relation 
S = (i)(N-l) 0. 

In Fig. 1(b) the v — i-FQHE gap is shown for n = 
and n — 1 LLs of the tachyonis system. For the n = LL, 
the FQHE gap does not depend on B and the gap exactly 
equals to the FQHE gap for the n = non-relativistic 
LL. This is because the n = tachyonic LL wave function 
consists of only the functions 4> m . For the n = 1 LL, the 
wave function is the B-dependent mixture of (j) m and 
4>i m . As a result the FQHE gap depends on the magnetic 
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FIG. 2: (a) The LL energy spectrum of a electron in a junction 
between two TIs is shown as a function of magnetic field for 
a few lowest LLs. The numbers next to the lines are the 
LL indexes. The TIs have the same Hamiltonian parameters, 
except A 1 , which takes the values: A 1 = 2.2 eV-A for TI-1 
and A 1 = 4.0 eV-A for TI-2. (b) The width in the z-direction 
of the LL wave functions is shown for n = and n = 1 LLs. 
The corresponding LL energy spectrum is shown in panel (a). 



field and changes from the n — non-relativistic LL value 
for B — > (point A) to the n = 1 graphene LL value for 
B = B* (point G) and finally to n = 1 non-relativistic 
LL value at point B (in the thermodynamic limit such a 
state becomes compressible). 

The maximum FQHE gap in a non-relativistic system 
corresponds to the green line in Fig. 1(b), while the max- 
imum FQHE gap in the graphene system corresponds to 
point G in Fig. 1(b). Therefore, comparing the data in 
Fig. 1(b), we conclude that within some range of the 
magnetic fields (which is shown in Fig. 1(b) by an oval 
curve), the FQHE gap in the model tachyonic system is 
the largest compared to other available 2D systems. The 
tachyonic dispersion relation provides an unique possi- 
bility to study, within a single tachyonic n = 1 LL, the 
properties of the n = non-relativistic LL [point A in 
Fig. 1(a)], n = 1 graphene LL [point G in Fig. 1(a)], and 
n = 1 non-relativistic LL [point B in Fig. 1(a)]. 

Three-dimensional (3D) model for the junction states 
between two TIs: Until now, we have discussed the 
magnetic field effects via an effective Hamiltonian for 
tachyons. The tachyonic dispersion can be realized in the 
junction of two TIs Q. The junction dispersion relation 
in this case is approximately described by the effective 
Hamiltonian ([1} . The realization of the tachyonic excita- 
tions as the junction states bring additional factors into 



consideration. For example, the junction states have a 
finite width in the z-direction which can reduce the in- 
teraction strength in that system. 

We consider the junction between two TI insulators: 
TI-1 for z < and TI-2 for z > 0. Here z = corresponds 
to the junction surface. The electronic properties of both 
TIs are described by the same type of low-energy effective 
3D Hamiltonian [9, [Q1 of the matrix form 



(A 2 /h)p + a x 



(A 2 /K)p_a x 
M(p)a z +iA 1 a x d z 



where d z — d/dz, and 

e(k) =C X - D x d 2 z 
M, 



(8) 



(9) 

M(k) = M + B.di - {B 2 /tf){pi+p 2 y ). (10) 
We assume that for TI-1 the constants in the above 



Hamiltonian are the same as for Bi 2 Se 3 [14j . while the 
for TI-2 the constants are different but close to the val- 
ues for Bi 2 Se 3 . We assume that only the constant A Y is 
different for the two TIs, i.e., A x = 2.2 eV-A for TI-1 
and A x = 4.0 eV-A for TI-2. All other constants in the 
Hamiltonian (JSJ) are kept the same for both TIs [4{. For 
these parameters, the junction states exhibit the tachy- 
onic dispersion [4|. The four-component wave function 
corresponding to the Hamiltonian ([5]) determines the am- 
plitudes of the wave functions at the positions of Bi and 
Se atoms: (Bi^, Se^, Bi^, Se^), where the arrows indicate 
the direction of the electron spin. 

The Hamiltonian of the TI in an external magnetic 
field, pointing along the z-direction, can be obtained 
from the Hamiltonian ([5]) by replacing the 2D momentum 
(p x ,p y ) by the generalized momentum (TT x ,n y ) [l5[ and 
introducing the Zeeman energy, A z = \g s [i B B. For the 
Hamiltonian ([5]) in a magnetic field, the wave function in 
the n-th LL has the general form 



X 2 {z)<t>n-l,m 



\ 
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Therefore, just as for the tachyonic states the wave func- 
tion ^„>{ is the mixture of n and n — 1 non-relativistic 



(n 

LL functions. For the n — LL, only % 3 
are non-zero. 

To find the LL junction states we follow the same 
procedure as for the LL surface states of a TI [Toh - 
\u\ . For each TI we find the general bulk solution of 

the Schrodinger equation in the form of ^ cx e A< >z , 
where A 1 -™ 1 ' is a complex constant, and m = 1 and 2 
for TI-1 and TI-2, respectively. This solution has a 
given energy, E, and a given LL index, n. The corre- 
sponding A*-" 1 -* are determined from a secular equation, 

det 



0) and xt 0) 



= 0, for each TI. For each en- 
ergy E and the LL index n, the secular equation pro- 
vides eight values of A,- (n, E), j = 1,...,8 with the 
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FIG. 3: The u=\ FQHE gap for n = (green line) and n = 1 
(red line) LLs in the junction of two TIs. The corresponding 
LLs are shown in Fig. 2(a) by green and red lines, respectively. 
The FQHE gap is calculated for a finite-size system with N = 
9 particles in spherical geometry with parameter S = 12. The 
FQHE gap is measured in units of Coulomb energy, e c = 



corresponding wave functions. Second, since we are look- 
ing for the localized LL junction states, we need to choose 
(for each TI) only four values of A,- out of eight with the 



properties: ReAj m) > for TI-1 (z < 0) and ReA}"* ; < 
for TI-2 (z > 0). We then choose the corresponding four 
wave functions (for each TI) as the basis and expand the 
solution for the LL junction state in this basis. Finally, 
the energy of the LL junction state is found from the con- 
dition of continuity of the wave function, \P„ (z), and 

( TT"t 

current [51-L Tl / 5k z \$>n in the junction. 

The spectrum of the LL junction states corresponding 
to the tachyonic dispersion is shown in Fig. 2(a). The 
spectrum is qualitatively similar to that [see Fig. 1(a)] 
obtained from the model 2D Hamiltonian. Both spec- 
tra have the finite range of magnetic fields and energies, 
where the LLs can be observed. At weak magnetic fields, 
the difference between the LL spectrum of the junction 
states and the 2D model is clearly visible. For a given 
LL index n, there are no junction states for weak mag- 
netic fields. These junction states are delocalized in the 
z-direction. To illustrate this delocalization we show in 
Fig. 2(b) the width of the n = and n = 1 LL wave 
functions in the z-direction. At a singular point of the LL 
spectrum, i.e., for B = B* , where the derivative of the LL 
energy with the magnetic field becomes infinitely large, 
the LL wave functions have the smallest width. This 
width increases with decreasing magnetic field and finally 
the LL junction states are delocalized in the z-direction. 
A similar behavior is observed for n = LL, but there are 



(m) 



no singular point in this case. Therefore, the LL energy 
spectrum of the junction states in the regime of tachy- 
onic dispersion can be well described within 2D effective 
model near the singular point B ~ B* . 

We have evaluated the FQHE gaps for the n = and 
n = 1 junction LLs. We have used the same approach as 
for the 2D model of the tachyonic excitations, discussed 
above. The results are shown in Fig. 3. Quantitatively 
the behavior of the FQHE gap as a function of the mag- 
netic field is similar to the 2D model of the tachyonic 
excitations [Fig. 1(b)]. Due to a finite width of the LL 
wave functions in the TI junction, there is a reduction of 
the inter-electron interaction strength and correspond- 
ingly the FQHE gap. This reduction is visible for n = 
LL, where a smaller FQHE gap and the magnetic field 
dependence of the FQHE gap is shown in Fig. 3. 

Although for the 2D tachyonic model the FQHE gap 
is the largest for the n — 1 LL, for the junction LLs the 
FQHE gap is the largest for the n — LL, due to the non- 
zero spin polarization of the n = 1 junction LL. This spin 
polarization is clearly visible from the general property 
of the LL wave function [Eq. (fTTj)]; only the components 

Xg™ -0 '' and xi" -0 "* of Vp£2o are non-zero and these com- 
ponents correspond to the spin-down polarization. The 
numerically found n = 1 LL wave functions also show 
partial spin-down polarization. As a result, the LL wave 
function have larger contribution from the 4> n =i non_ 
relativistic LL function, which reduces the inter-electron 
interaction strength and the FQHE gap. 

To summarize: we have investigated the magnetic field 
effects of tachyonic excitations along the interface of 
two topological insulators. We used an effective two- 
dimensional model Hamiltonian for tachyons and the 
three-dimensional model for the junction states of the 
two TIs, both developed by us The Landau levels in 
both these models show very similar behaviors. Unlike 
in graphene or in conventional electron systems, only a 
few LLs are found to exist for the tachyons. Only one 
LL (n = 0) survives for large magnetic fields. The v = | 
FQHE is the strongest (within a limited range of the mag- 
netic field) when compared with that for conventional 
electron systems and graphene. Interestingly, the FQHE 
in the n = 1 LLs for tachyons describes the FQHE of the 
n = 0, 1 LLs of the non-relativistic electron system and 
that of the n — 1 graphene LL in different regions of the 
magnetic field. Experimental confirmation of these prop- 
erties of the Landau levels would provide strong evidence 
on the existence of elusive tachyons. 
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